The Na(p, 7) Mg Reaction and Oxygen-Neon Novae 
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The 21 Na(p, 7) 22 Mg reaction is expected to play an important role in the nucleosynthesis of 22 Na 
in Oxygen-Neon novae. The decay of 22 Na leads to the emission of a characteristic 1.275 MeV 
gamma-ray line. This report provides the first direct measurement of the rate of this reaction 
using a radioactive 21 Na beam, and discusses its astrophysical implications. The energy of the 
important state was measured to be E c . m .= 205.7 ± 0.5 keV with a resonance strength ui-y — 
1.03 ± 0.16 s tat ± 0.14 SHS meV. 

PACS numbers: 25.60.-t, 26.30.+k 
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The synthesis of light and intermediate-mass elements 
can take place through radiative proton captures on un- 
stable nuclei during explosive stellar events. One as- 
trophysical site where such processes can occur involves 
classical novae, stellar explosions powered by thermonu- 
clear runaways on accreting ONe or CO white dwarf 
stars 0, 0, E|- I n this hydrogen burning process, nu- 
clear activity involves different cycles, depending on the 
nova type and on the temperatures achieved during the 
explosion. A predominant nuclear activity in ONe no- 
vae takes place in the NeNa cycle, initiated by radiative 
proton captures on the abundant seed nuclei 20 Ne. 

Nucleosynthesis in the NeNa cycle during nova out- 
bursts leads to the synthesis of the astronomically im- 
portant, but unstable 22 Na nucleus. Its /3-decay (ti/2 
— 2.6 yr) leads to the emission of a 1.275 MeV 7-ray, 
following population of the first excited state of 22 Ne. 
In fact, this 7-ray is an ideal observable for nova events 
as first suggested by Clayton & Hoyle Q. Thus far, 
observational searches performed with NASA's COMP- 
TEL on-board CGRO satellite of five ONe novae have 
not found this 7-ray signature 0|. Whereas the inferred 
upper limits are in agreement with recent results from 
ONe nova models 0, L| , the reduction of the nuclear un- 
certainties associated with the main reactions involved in 
the synthesis of 22 Na is critically important in order to 
predict how much 22 Na can be produced in a typical nova 



event, and at what distance a nova explosion is expected 
to provide a detectable flux of 7-rays. 

Another aspect that stresses the astronomical inter- 
est of 22 Na relies on the identification of presolar grains 
likely condensed in the ejecta from nova outbursts. Tra- 
ditionally, they have been identified by low 20 Ne/ 22 Ne 
ratios (where 22 Ne is attributed to in-situ 22 Na decay). A 
22 Na/C ratio of 9 x 10 -6 |6j has been determined recently 
in the graphite grain KFBla-161, in which other isotopic 
ratios resemble those found in the envelopes ejected by 
nova outbursts. Again, a more accurate knowledge of re- 
actions in the synthesis of 22 Na in novae would further 
assist in identifying presolar grains from novae and for 
tuning models accordingly. 

Synthesis of 22 Na in novae takes place following two 
possible reaction paths (Fig.GJ: in the first ("cold" NeNa 
cycle), 21 Na forms from the seed 20 Ne which then leads to 
21 Na(/3 + ) 21 Ne(p,7) 22 Na; in the second path, associated 
with higher temperatures ("hot" NeNa cycle), proton- 
capture on 21 Na dominates over its /3-decay, followed by 
21 Na(p,7) 22 Mg(/3+) 22 Na. There is little net mass flow 
from 22 Mg to 23 Al due to the low Q-value for photo- 
disintegration of 23 Al 0- Current models of ONe no- 
vae indicate that the unknown rate of 21 Na(p,7) is the 
main source of uncertainty associated with calculating 
the amount of 22 Na in nova outbursts [a, 111- The purpose 
of this paper is to report on the first direct measurement 
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FIG. 1: The combined "cold" and "hot" NeNa reaction cycles. 
The isotope 21 Na will either /3-decay into 21 Ne(the "cold" 
NeNa cycle) or capture a proton leading to 22 Mg (the "hot" 
cycle) depending upon the temperature and the reaction rate. 



of this rate. 

Under nova conditions the capture reaction rate, 
Na(&v), is expected to be dominated by one or more 
narrow resonances. Each resonance contributes to the re- 
action rate in direct proportion to its resonance strength, 
wy, and depends exponentially on the resonance energy, 



Er. In units of cnrs mol , it is given by, 



N A (av) = 1.54 x 10 n (/iT9)- 3/2 
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with Na Avogadro's number, (i the reduced mass in u, Tg 
the temperature in units of GK, (av) the thermally aver- 
aged nuclear cross section, and u>j and Er in MeV fici| . 
The narrow resonance thick target yield, Y , at maximum 

is EH, 
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with A the centre-of-mass de Broglie wave length, M the 
(heavy) projectile nucleus mass, m the (light) target nu- 
cleus mass, and 4j| the energy loss per atom/cm 2 (lab). 
Thus, measurement of the maximum thick target yield 
can determine the resonance strength, wy. 

Figure shows the 22 Mg level scheme Uj,H^. Calcu- 
lation of the Gamow window indicates that the 212 keV, 
t = resonance will be the dominant contributor (as 
compared to other higher resonances and direct capture) 
to the 21 Na(p, 7) 22 Mg reaction at all nova temperatures 
from 0.2 to 0.35 GK. We report here a measurement of 
the strength, wy, and a revised energy, Er, for this res- 
onance in the 21 Na(p, 7) 22 Mg reaction. 

The experiment was carried out at the TRIUMF-ISAC 
radioactive beams facility located in Vancouver, Canada. 
Fifteen /iA of 500 MeV protons bombarded a thick target 
of SiC resulting in an intense (~ 10 9 s _1 ), pure (~ 100%) 
21 Na beam extracted from a surface ion source and a 
high resolution mass analyzer 0] . It was accelerated us- 
ing the new ISAC linear accelerator, resulting in beams 
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FIG. 2: The Mg level scheme of those states of astrophysical 
interest for ONe nova, shown with solid lines [l2T UsT ] . The 
numbers on the far left denote centre-of-mass energies (E x - Q) 
in units of keV. The state at 5.837 MeV was observed once 
but not confirmed in other studies \VA H3L . 



with energies variable from 0.15 to 1.5 MeV/u The 
study was performed using inverse kinematics with the 
DRAGON (Detector of Recoils And Gammas Of Nuclear 
reactions) facility. DRAGON consists of a windowless 
gas target (effective length of 12.3 cm) surrounded by a 
gamma array (30 units of BGO), and followed by a two- 
stage, recoil mass separator, 21 m in length (from tar- 
get centre to focal plane). Separation of the rare recoil 
from more intense beam is achieved using magnetic and 
electric dipoles. Following an initial selection of a single 
(optimal) charge state |ljfl in the first magnetic dipole, 
energy dispersion in the electric dipole allows mass sepa- 
ration, and the process is repeated in the second stage. A 
DSSSD (Double Sided Silicon Strip Detector) was used 
at the focal plane of DRAGON to detect the 22 Mg re- 
coils. A more complete description of DRAGON can be 
found elsewhere 0, 0] . 

A radioactive beam of 21 Na (q=5 + ) at typical inten- 
sities up to 5 x 10 8 s" 1 was delivered to the DRAGON 
hydrogen gas target (4.6 Torr). The gas target received 
a total of ~ 10 13 21 Na atoms for this study. Data taking 
was done in both singles and coincidence modes; the co- 
incidence mode required a "start" timing signal from the 
7-array in coincidence with a "stop" timing signal from 
the DSSSD. Figure EH shows resonant-capture spectra for 
a beam energy of 220 keV/u. Counts within the box in 
Fig. EK were considered to be valid capture events. Their 
recoil energy distribution is presented in Fig. Eh. Fig. 0t 
is the recoil time-of-fiight spectrum for events satisfying 
the cut on gamma-ray energy. The distribution of the 
hit BGO detector position along the beam axis (Fig. 01) 
shows that the resonance was near the centre of the gas 
target at beam energy 220 keV/u {E c .m.= 211 keV). 

The beam energies were measured by adjusting the 
field of the first magnetic dipole in the separator so as to 
position the beam on the ion-optical axis at an energy- 
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FIG. 3: Resonant-capture spectra for a 21 Na beam energy of 
220 keV/u. (a) Valid events enclosed by a two-dimensional- 
cut box above a background of random-coincidence events, (b) 
the recoil-energy distribution of the events in the DSSSD, 
selected by the box in (a), (c) the recoil TOF distribution for 
events above the 7-ray threshold energy, (d) distribution of 
box-selected 7-ray events observed in the BGO array along 
the target length, with a Gaussian fit. 



dispersed focus. Using the design bending radius of the 
dipole (1 m), it was possible to calculate beam energy in 
terms of the dipole field. The expected relationship was 
confirmed by measuring a number of known resonances 
with stable beams. The lower panel of Fig. 0] shows the 
yield curve for one of these studies, the 24 Mg(p, 7) 25 A1 
reaction, demonstrating our agreement (inflection point 
of 214.4 ± 0.5 keV) with the literature resonance energy 
of 214.0 ± 0.1 keV 20]. As shown in the upper panel, we 
find the energy (inflection point) for the 21 Na(p,7) 22 Mg 
resonance to be 205.7 ± 0.5, and not 212 keV (see Fig.[2J), 
the difference between the Q value and the level exci- 
tation energy, 5713.9 ± 1.2 keV Given that the 
latter value is based upon a direct gamma de-excitation 
measurement of the 5713.9 keV level, this disagreement 
could be explained by a modified mass excess for 22 Mg; 
our data imply a value of -403.2 ± 1.3 keV rather than 
-396.8 keV Q. 

Figure 0] (upper panel) shows the thick target yield 
curve corrected/scaled for various factors listed in Ta- 
ble HJ The efficiency of the BGO array as a function of 
7-ray energy and resonance position in the targ et was 
calculated using the GEANT program HI Hj. The 
variation of resonance position with beam energy re- 
sulted in the following calculated efficiencies: 45% for 
202 keV < E < 207 keV, 48% at 211 keV, and 46% 
above 216 keV. The systematic error was deduced from 
values of the array efficiency measured with stable beam 
reactions. The separator transmission (98%) and 
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FIG. 4: The upper panel displays the thick target yield data 
for the 21 Na(p, 7) 22 Mg reaction, with the solid line show- 
ing the nominal target thickness for 4.6 Torr. Yield of the 
24 Mg(p,7) 25 Al reaction for the resonance at E c . m .=214 keV, 
used for beam energy calibration, is displayed in the lower 
panel. Statistical errors only are displayed in both. 



TABLE I: Summary of systematic errors. 



Factors 



Value Syst.Error(%) 



BGO array efficiency (Q211 keV) 
Separator transmission 

DSSSD efficiency 
Charge state fraction 
Integrated beam (@211 keV) 
dE/dx (eV/(atom/cm 2 ))i ai) 



0.48 
0.98 
0.99 
0.44 
3.62 x 10 13 
8.18 x 10~ 14 
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1 

3 
4 
5 



DSSSD detection efficiency (99%) [25| were determined 
separately, and the fraction of the charge state selected 
(44%) was measured with a 24 Mg beam of 220 keV/u. 
At 4.6 Torr, charge state equilibrium in H2 gas was mea- 
sured to be attained within 2 mm |l7j . The energy loss 
in the target (4.6 Torr) was measured to be 14.4 keV/u 
(lab) or 8.18 x 10 _14 eV/(atom/cm 2 ), in agreement with 

SRIM m. 

The data of Fig. 0] (upper panel) were obtained by 
maximum likelihood combination of several runs at each 
energy [27j . The error bars on the zero counts seen at 
off- resonance energies are 68% confidence limits. Table H] 
presents a summary of systematic errors. Using Eq. 
and only the mid-target data point (211 keV), a yield 
of (5.76 ± 0.88) x 10" 12 per incident 21 Na, results in a 
resonance strength of W7 = 1.03±0.16 s t a t±0.14 sys meV. 

The effect of these results on the calculated stellar re- 
action rate is shown in Fig. The rate is reduced over 
that determined by shell model calculations of wyas re- 
ported in 01 1 an d enhanced over that found in ||. An 
analysis of the impact of the new measurements on the 
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FIG. 5: The stellar rate for the 21 Na(p, 7) 22 Mg reaction using 
Eq. 2 with typical novae temperatures and our measured 
values for 017 and Eb.— 0.206 MeV (solid line with hatched 
area reflecting errors), in comparison with other works; upper 
curve |T2] and lower curve ||. 



synthesis of 22 Na in novae was performed. A new model 
of a nova outburst, using an ONe white dwarf of 1.25 
solar mass, has been computed from the onset of accre- 
tion up to the explosion and ejection stages, by means 
of a spherically symmetric, implicit, hydrodynamic code, 
in Lagrangian formulation (see |3j for details). Results 
have been compared with a model evolved with the pre- 
vious prescription of the 21 Na(p,7) 22 Mg rate §|. As a 
result of the higher contribution of the 5.714 MeV level 
(Fig. [5}, a slightly lower amount of 22 Na (a mean mass 
fraction of 2.8 x 1CP 4 , compared with the previous esti- 
mate of 3.5 x 1CP 4 ) is found. The small decrease in the 
22 Na yield results from the fact that increasing the proton 
capture rate on 21 Na favours the synthesis path through 
21 Na(p,7) 22 Mg(/3 + ) 22 Na, hence reducing the role of the 
alternative 21 Na(/3 + ) 21 Ne(p,7) 22 Na path. In these newly 
derived conditions of increased proton capture on 21 Na, 
22 Na production takes place earlier in the outburst, at 
a time when the envelope has not yet significantly ex- 
panded and cooled down (contrary to the case when a 
lower 21 Na(p,7) rate is adopted), and hence the temper- 
ature in the envelope is still high enough to allow proton 
captures on 22 Na, that reduce its final content in the 
ejecta. 

Up to now, 7-ray flux determinations were limited by 
a large uncertainty in the the 21 Na(p,7) and 22 Na(p,7) 
rates, which translated into an overall uncertainty in 
the 22 Na yields of a factor of ~ 3. The maximum 
detectability distance was, accordingly, uncertain by a 
factor of ~ 2. Such uncertainty, mainly due to the 
previously unknown reaction rate, has been largely re- 
duced with the present experimental determination of 
CJ7 = 1.034=0. 16 s tat±0.14 S j /s meV. These results provide 



a firmer basis for predictions of the expected gamma-ray 
signature at 1.275 MeV associated with 22 Na decay in 
ONe novae, and confirm the previous determination of 
1 kiloparsec for a typical ONe nova [2^, [2j| observed 
with ESA's (European Space Agency) INTEGRAL spec- 
trometer, SPI. Furthermore, the smaller uncertainty in 
the rate also indicates that the predicted 22 Na yields are 
not in conflict with the upper limits derived from several 
observational searches. 
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